Extended-volume image-derived models of coronary microcirculation.

OBJECTIVE: Recent advances in tissue clearing and high-throughput imaging have enabled the acquisition of extended-volume microvasculature images at a submicron resolution. The objective of this study was to extract information from this type of images by integrating a sequence of 3D image processing steps on Terabyte scale datasets. METHODS: We acquired coronary microvasculature images throughout an entire short-axis slice of a 3-month-old Wistar-Kyoto rat heart. This dataset covered 13 × 10 × 0.6 mm at a resolution of 0.933 × 0.933 × 1.866 µm and occupied 700 Gigabytes of disk space. We used chunk-based image segmentation, combined with an efficient graph generation technique, to quantify the microvasculature in the large-scale images. Specifically, we focused on the microvasculature with a vessel diameter up to 15 µm. RESULTS: Morphological data for the complete short-axis ring were extracted within 16 h using this pipeline. From the analyses, we identified that microvessel lengths in the rat coronary microvasculature varied from 6 to 300 µm. However, their distribution was heavily skewed toward shorter lengths, with a mode of 16.5 µm. In contrast, vessel diameters ranged from 3 to 15 µm and had an approximately normal distribution of 6.5 ± 2 µm. CONCLUSION: The tools and techniques from this study will serve other investigations into the microcirculation, and the wealth of data from this study will enable the analysis of biophysical mechanisms using computer models.

Personalized surgical planning for coronary bypass graft configurations using patient-specific computational modeling to avoid flow competition in arterial grafts.

Objectives: Flow competition between coronary artery bypass grafts (CABG) and native coronary arteries is a significant problem affecting arterial graft patency. The objectives of this study were to compare the predictive hemodynamic flow resulting from various total arterial grafting configurations and to evaluate whether the use of computational fluid dynamics (CFD) models capable of predicting flow can assist surgeons to make better decisions for individual patients by avoiding poorly functioning grafts. Methods: Sixteen cardiac surgeons declared their preferred CABG configuration using bilateral internal mammary and radial arteries for each of 5 patients who had differing degrees of severe triple vessel coronary disease. Surgeons selected both a preferred 'aortic' strategy, with at least one graft arising from the ascending aorta, and a preferred "anaortic" strategy which could be performed as a "no-aortic touch" operation. CT coronary angiograms of the 5 patients were coupled to CFD models using a novel flow solver "COMCAB." Twelve different CABG configurations were compared for each patient of which 4 were "aortic" and 8 were "anaortic." Surgeons then selected their preferred grafting configurations after being shown predictive hemodynamic metrics including functional assessment of stenoses (instantaneous wave-free ratio; fractional flow reserve), transit time flowmetry graft parameters (mean graft flow; pulsatility index) and myocardial perfusion. Results: A total of 87.5% (7/8) of "anaortic" configurations compared to 25% (1/4) of "aortic" configurations led to unsatisfactory grafts in at least 1 of the 5 patients (P = 0.038). The use of the computational models led to a significant decrease in the selection of unsatisfactory grafting configurations when surgeons employed "anaortic" (21.25% (17/80) vs. 1.25% (1/80), P < 0.001) but not "aortic" techniques (5% (4/80) vs. 0% (0/80), P = 0.64). Similarly, there was an increase in the selection of ideal configurations for "anaortic" (6.25% (5/80) vs. 28.75% (23/80), P < 0.001) but not "aortic" techniques (65% (52/80) vs. 61.25% (49/80), P = 0.74). Furthermore, surgeons who planned to use more than one unique "anaortic" configuration across all 5 patients increased (12.5% (2/16) vs. 87.5% (14/16), P<0.001). Conclusions: "COMCAB" is a promising tool to improve personalized surgical planning particularly for CABG configurations involving composite or sequential grafts which are used more frequently in anaortic operations.

A predictive patient-specific computational model of coronary artery bypass grafts for potential use by cardiac surgeons to guide selection of graft configurations.

Cardiac surgeons face a significant degree of uncertainty when deciding upon coronary artery bypass graft configurations for patients with coronary artery disease. This leads to significant variation in preferred configuration between different surgeons for a particular patient. Additionally, for the majority of cases, there is no consensus regarding the optimal grafting strategy. This situation results in the tendency for individual surgeons to opt for a "one size fits all" approach and use the same grafting configuration for the majority of their patients neglecting the patient-specific nature of the diseased coronary circulation. Quantitative metrics to assess the adequacy of coronary bypass graft flows have recently been advocated for routine intraoperative use by cardiac surgeons. In this work, a novel patient-specific 1D-0D computational model called "COMCAB" is developed to provide the predictive haemodynamic parameters of functional graft performance that can aid surgeons to avoid configurations with grafts that have poor flow and thus poor patency. This model has significant potential for future expanded applications.

Venous thrombosis risk during and after medical and surgical hospitalizations: The medical inpatient thrombosis and hemostasis (MITH) study.

BACKGROUND: Thirty to seventy percent of all venous thromboembolism (VTE) events are associated with hospitalization. The absolute and relative risks during and after hospitalization are poorly characterized. OBJECTIVES: Quantify the absolute rate and relative risk of VTE during and up to 3 months after medical and surgical hospitalizations. PATIENTS/METHODS: We conducted an observational cohort study between 2010 and 2016 of patients cared for by the University of Vermont (UVM) Health Network's primary care population. Cox proportional hazard models with hospitalization modeled as a time-varying covariate were used to estimate VTE risk. RESULTS: Over 4.3 years of follow-up, 55 220 hospitalizations (156 per 1000 person-years) and 713 first venous thromboembolism events (2.0 per 1000 person-years) occurred. Among individuals not recently hospitalized, the rate of venous thromboembolism was 1.4 per 1000 person-years and 71.8 per 1000 person-years during hospitalization. During the first, second, and third months after discharge, the rates of venous thromboembolism were 35.1, 11.3, and 5.2 per 1000 person-years, respectively. Relative to those not recently hospitalized, the age- and sex-adjusted HRs of venous thromboembolism were 38.0 (95% CI 28.0, 51.5) during hospitalization, and 18.4 (95% CI 15.0, 22.6), 6.3 (95% CI 4.3, 9.0), and 3.0 (95% CI 1.7, 5.4) during the first, second, and third months after discharge, respectively. Stratified by medical versus surgical services the rates were similar. CONCLUSION: Hospitalization and up to 3 months after discharge were strongly associated with increased venous thromboembolism risk. These data quantify this risk for use in future studies.

It's clearly the heart! Optical transparency, cardiac tissue imaging, and computer modelling.

Recent developments in clearing and microscopy enable 3D imaging with cellular resolution up to the whole organ level. These methods have been used extensively in neurobiology, but their uptake in other fields has been much more limited. Application of this approach to the human heart and effective use of the data acquired present challenges of scale and complexity. Four interlinked issues need to be addressed: 1) efficient clearing and labelling of heart tissue, 2) fast microscopic imaging of human-scale samples, 3) handling and processing of multi-terabyte 3D images, and 4) extraction of structural information in computationally tractable structure-based models of cardiac function. Preliminary studies show that each of these requirements can be achieved with the appropriate application and development of existing technologies.

Effects of Thyroid Function on Hemostasis, Coagulation, and Fibrinolysis: A Mendelian Randomization Study.

Background: Untreated hypothyroidism is associated with acquired von Willebrand syndrome, and hyperthyroidism is associated with increased thrombosis risk. However, the causal effects of thyroid function on hemostasis, coagulation, and fibrinolysis are unknown. Methods: In a two-sample Mendelian randomization (MR) study with genome-wide association variants, we assessed causality of genetically predicted hypothyroidism (N = 134,641), normal-range thyrotropin (TSH; N = 54,288) and free thyroxine (fT4) (N = 49,269), hyperthyroidism (N = 51,823), and thyroid peroxidase antibody positivity (N = 25,821) on coagulation (activated partial thromboplastin time, von Willebrand factor [VWF], factor VIII [FVIII], prothrombin time, factor VII, fibrinogen) and fibrinolysis (D-dimer, tissue plasminogen activator [TPA], plasminogen activator inhibitor-1) from the CHARGE Hemostasis Consortium (N = 2583-120,246). Inverse-variance-weighted random effects were the main MR analysis followed by sensitivity analyses. Two-sided p < 0.05 was nominally significant, and p < 0.0011[ = 0.05/(5 exposures × 9 outcomes)] was Bonferroni significant for the main MR analysis. Results: Genetically increased TSH was associated with decreased VWF [ß(SE) = -0.020(0.006), p = 0.001] and with decreased fibrinogen [ß(SE) = -0.008(0.002), p = 0.001]. Genetically increased fT4 was associated with increased VWF [ß(SE) = 0.028(0.011), p = 0.012]. Genetically predicted hyperthyroidism was associated with increased VWF [ß(SE) = 0.012(0.004), p = 0.006] and increased FVIII [ß(SE) = 0.013(0.005), p = 0.007]. Genetically predicted hypothyroidism and hyperthyroidism were associated with decreased TPA [ß(SE) = -0.009(0.024), p = 0.024] and increased TPA [ß(SE) = 0.022(0.008), p = 0.008], respectively. MR sensitivity analyses showed similar direction but lower precision. Other coagulation and fibrinolytic factors were inconclusive. Conclusions: In the largest genetic studies currently available, genetically increased TSH and fT4 may be associated with decreased and increased synthesis of VWF, respectively. Since Bonferroni correction may be too conservative given the correlation between the analyzed traits, we cannot reject nominal associations of thyroid traits with coagulation or fibrinolytic factors.

Determining the optimal dosing of a novel combination regimen of ceftazidime/avibactam with aztreonam against NDM-1-producing Enterobacteriaceae using a hollow-fibre infection model.

BACKGROUND: MBL-producing strains of Enterobacteriaceae are a major public health concern. We sought to define optimal combination regimens of ceftazidime/avibactam with aztreonam in a hollow-fibre infection model (HFIM) of MBL-producing strains of Escherichia coli and Klebsiella pneumoniae. METHODS: E. coli ARLG-1013 (blaNDM-1, blaCTX-M, blaCMY, blaTEM) and K. pneumoniae ARLG-1002 (blaNDM-1, blaCTXM-15, blaDHA, blaSHV, blaTEM) were studied in the HFIM using simulated human dosing regimens of ceftazidime/avibactam and aztreonam. Experiments were designed to evaluate the effect of staggered versus simultaneous administration, infusion duration and aztreonam daily dose (6 g/day versus 8 g/day) on bacterial killing and resistance suppression. Prospective validation experiments for the most active combination regimens were performed in triplicate to ensure reproducibility. RESULTS: Staggered administration of the combination (ceftazidime/avibactam followed by aztreonam) was found to be inferior to simultaneous administration. Longer infusion durations (2 h and continuous infusion) also resulted in enhanced bacterial killing relative to 30 min infusions. The rate of killing was more pronounced with 8 g/day versus 6 g/day aztreonam combination regimens for both tested strains. In the prospective validation experiments, ceftazidime/avibactam with aztreonam dosed every 8 and 6 h, respectively (ceftazidime/avibactam 2/0.5 g every 8 h + aztreonam 2 g every 6 h), or ceftazidime/avibactam with aztreonam as continuous infusions resulted in maximal bacterial killing and resistance suppression over 7 days. CONCLUSIONS: Simultaneous administration of aztreonam 8 g/day given as a continuous or 2 h infusion with ceftazidime/avibactam resulted in complete bacterial eradication and resistance suppression. Further study of this combination is needed with additional MBL-producing Gram-negative pathogens. The safety of this double ß-lactam strategy also warrants further study in Phase 1 clinical trials.

Reconstruction of coronary circulation networks: A review of methods.

Building anatomically accurate models of the coronary vascular system enables potentially deeper understandings of coronary circulation. To achieve this, (a) images at different levels of vascular network-arteries, arterioles, capillaries, venules, and veins-need to be obtained through suitable imaging modalities; and (b) from images, morphological and topological information needs to be extracted using image processing techniques. While there are several modalities that enable the imaging of large vessels, microcirculation imaging-capturing vessels having diameter lesser than 100 µm-has to date been typically confined to small regions of the heart. This spatially limited microcirculatory information has often been used within cardiac models, with the potentially erroneous assumption that it is representative of the whole organ. However, with the recent advancements in imaging and image processing, it is rapidly becoming feasible to acquire, process, and quantify microcirculation data at the scale of whole organ. In this review, we summarize the progress toward this goal followed through a presentation of the current state-of-the-art imaging and image processing techniques in the context of coronary microcirculation extraction, prominently but not exclusively, from small animals.

Electrically reconfigurable terahertz signal processing devices using liquid metal components.

Many applications of terahertz (THz) technology require the ability to actively manipulate a free space THz beam. Yet, although there have been many reports on the development of devices for THz signal processing, few of these include the possibility of electrical control of the functionality, and novel ideas are needed for active and reconfigurable THz devices. Here, we introduce a new approach, based on the integration of electrically actuated liquid metal components in THz waveguides. This versatile platform offers many possibilities for control of THz spectral content, wave fron"ts, polarization, and power flow. We demonstrate two illustrative examples: the first active power-splitting switch, and the first channel add-drop filter. We show that both of these devices can be used to electrically switch THz communication signals while preserving the information in a high bit-rate-modulated data stream.

Decreasing Compensatory Ability of Concentric Ventricular Hypertrophy in Aortic-Banded Rat Hearts.

The cardiac system compensates for variations in physiological and pathophysiological conditions through a dynamic remodeling at the organ, tissue, and intracellular levels in order to maintain function. However, on longer time scales following the onset of ventricular pressure overload, such remodeling may begin to inhibit physiological function and ultimately lead to heart failure. This progression from compensatory to decompensatory behavior is poorly understood, in particular owing to the absence of a unified perspective of the concomitantly remodeling subsystems. To address this issue, the present study investigates the evolution of compensatory mechanisms, in response to overload, by integrating diffusion-tensor MRI, echocardiography, and intracellular and hemodynamic measurements within consistent computational simulations of aortic-banded rat hearts. This approach allows a comparison of the relative leverage of different cardiac properties (geometry, passive mechanical stiffness, fiber configuration, diastolic and peak calcium concentrations, calcium-binding affinity, and aortic impedance) to affect cardiac contraction. Measurements indicate that, following aortic banding, an ejection fraction (EF) of 75% was maintained, relative to control rats, despite significant remodeling of the left-ventricular wall thickness (increasing by ~90% over 4 weeks). Applying our framework, we identified the left-ventricular wall thickness (concentric hypertrophy) and the intracellular calcium dynamics as playing the dominant roles in preserving EF acutely, whereas the significance of hypertrophy decreased subsequently. This trend suggests an increasing reliance on intracellular mechanisms (average increase ~50%), rather than on anatomical features (average decrease ~60%), to achieve compensation of pump function in the early phase of heart failure.

Compensatory and decompensatory alterations in cardiomyocyte Ca2+ dynamics in hearts with diastolic dysfunction following aortic banding.

KEY POINTS: At the cellular level cardiac hypertrophy causes remodelling, leading to changes in ionic channel, pump and exchanger densities and kinetics. Previous studies have focused on quantifying changes in channels, pumps and exchangers without quantitatively linking these changes with emergent cellular scale functionality. Two biophysical cardiac cell models were created, parameterized and validated and are able to simulate electrophysiology and calcium dynamics in myocytes from control sham operated rats and aortic-banded rats exhibiting diastolic dysfunction. The contribution of each ionic pathway to the calcium kinetics was calculated, identifying the L-type Ca2+ channel and sarco/endoplasmic reticulum Ca2+ ATPase as the principal regulators of systolic and diastolic Ca2+ , respectively. Results show that the ability to dynamically change systolic Ca2+ , through changes in expression of key Ca2+ modelling protein densities, is drastically reduced following the aortic banding procedure; however the cells are able to compensate Ca2+ homeostasis in an efficient way to minimize systolic dysfunction. ABSTRACT: Elevated left ventricular afterload leads to myocardial hypertrophy, diastolic dysfunction, cellular remodelling and compromised calcium dynamics. At the cellular scale this remodelling of the ionic channels, pumps and exchangers gives rise to changes in the Ca2+ transient. However, the relative roles of the underlying subcellular processes and the positive or negative impact of each remodelling mechanism are not fully understood. Biophysical cardiac cell models were created to simulate electrophysiology and calcium dynamics in myocytes from control rats (SHAM) and aortic-banded rats exhibiting diastolic dysfunction. The model parameters and framework were validated and the fitted parameters demonstrated to be unique for explaining our experimental data. The contribution of each ionic pathway to the calcium kinetics was calculated, identifying the L-type Ca2+ channel (LCC) and the sarco/endoplasmic reticulum Ca2+ -ATPase (SERCA) as the principal regulators of systolic and diastolic Ca2+ , respectively. In the aortic banding model, the sensitivity of systolic Ca2+ to LCC density and diastolic Ca2+ to SERCA density decreased by 16-fold and increased by 23%, respectively, relative to the SHAM model. The energy cost of ionic homeostasis is maintained across the two models. The models predict that changes in ionic pathway densities in compensated aortic banding rats maintain Ca2+ function and efficiency. The ability to dynamically alter systolic function is significantly diminished, while the capacity to maintain diastolic Ca2+ is moderately increased.

A model of cardiac contraction based on novel measurements of tension development in human cardiomyocytes.

Experimental data from human cardiac myocytes at body temperature is crucial for a quantitative understanding of clinically relevant cardiac function and development of whole-organ computational models. However, such experimental data is currently very limited. Specifically, important measurements to characterize changes in tension development in human cardiomyocytes that occur with perturbations in cell length are not available. To address this deficiency, in this study we present an experimental data set collected from skinned human cardiac myocytes, including the passive and viscoelastic properties of isolated myocytes, the steady-state force calcium relationship at different sarcomere lengths, and changes in tension following a rapid increase or decrease in length, and after constant velocity shortening. This data set is, to our knowledge, the first characterization of length and velocity-dependence of tension generation in human skinned cardiac myocytes at body temperature. We use this data to develop a computational model of contraction and passive viscoelasticity in human myocytes. Our model includes troponin C kinetics, tropomyosin kinetics, a three-state crossbridge model that accounts for the distortion of crossbridges, and the cellular viscoelastic response. Each component is parametrized using our experimental data collected in human cardiomyocytes at body temperature. Furthermore we are able to confirm that properties of length-dependent activation at 37°C are similar to other species, with a shift in calcium sensitivity and increase in maximum tension. We revise our model of tension generation in the skinned isolated myocyte to replicate reported tension traces generated in intact muscle during isometric tension, to provide a model of human tension generation for multi-scale simulations. This process requires changes to calcium sensitivity, cooperativity, and crossbridge transition rates. We apply this model within multi-scale simulations of biventricular cardiac function and further refine the parametrization within the whole organ context, based on obtaining a healthy ejection fraction. This process reveals that crossbridge cycling rates differ between skinned myocytes and intact myocytes.

Accurate and Standardized Coronary Wave Intensity Analysis.

OBJECTIVE: Coronary wave intensity analysis (cWIA) has increasingly been applied in the clinical research setting to distinguish between the proximal and distal mechanical influences on coronary blood flow. Recently, a cWIA-derived clinical index demonstrated prognostic value in predicting functional recovery postmyocardial infarction. Nevertheless, the known operator dependence of the cWIA metrics currently hampers its routine application in clinical practice. Specifically, it was recently demonstrated that the cWIA metrics are highly dependent on the chosen Savitzky-Golay filter parameters used to smooth the acquired traces. Therefore, a novel method to make cWIA standardized and automatic was proposed and evaluated in vivo. METHODS: The novel approach combines an adaptive Savitzky-Golay filter with high-order central finite differencing after ensemble-averaging the acquired waveforms. Its accuracy was assessed using in vivo human data. The proposed approach was then modified to automatically perform beat wise cWIA. Finally, the feasibility (accuracy and robustness) of the method was evaluated. RESULTS: The automatic cWIA algorithm provided satisfactory accuracy under a wide range of noise scenarios (≤10% and ≤20% error in the estimation of wave areas and peaks, respectively). These results were confirmed when beat-by-beat cWIA was performed. CONCLUSION: An accurate, standardized, and automated cWIA was developed. Moreover, the feasibility of beat wise cWIA was demonstrated for the first time. SIGNIFICANCE: The proposed algorithm provides practitioners with a standardized technique that could broaden the application of cWIA in the clinical practice as enabling multicenter trials. Furthermore, the demonstrated potential of beatwise cWIA opens the possibility investigating the coronary physiology in real time.

Beyond Bernoulli: Improving the Accuracy and Precision of Noninvasive Estimation of Peak Pressure Drops.

BACKGROUND: Transvalvular peak pressure drops are routinely assessed noninvasively by echocardiography using the Bernoulli principle. However, the Bernoulli principle relies on several approximations that may not be appropriate, including that the majority of the pressure drop is because of the spatial acceleration of the blood flow, and the ejection jet is a single streamline (single peak velocity value). METHODS AND RESULTS: We assessed the accuracy of the Bernoulli principle to estimate the peak pressure drop at the aortic valve using 3-dimensional cardiovascular magnetic resonance flow data in 32 subjects. Reference pressure drops were computed from the flow field, accounting for the principles of physics (ie, the Navier-Stokes equations). Analysis of the pressure components confirmed that the spatial acceleration of the blood jet through the valve is most significant (accounting for 99% of the total drop in stenotic subjects). However, the Bernoulli formulation demonstrated a consistent overestimation of the transvalvular pressure (average of 54%, range 5%-136%) resulting from the use of a single peak velocity value, which neglects the velocity distribution across the aortic valve plane. This assumption was a source of uncontrolled variability. CONCLUSIONS: The application of the Bernoulli formulation results in a clinically significant overestimation of peak pressure drops because of approximation of blood flow as a single streamline. A corrected formulation that accounts for the cross-sectional profile of the blood flow is proposed and adapted to both cardiovascular magnetic resonance and echocardiographic data.

Non-invasive Model-Based Assessment of Passive Left-Ventricular Myocardial Stiffness in Healthy Subjects and in Patients with Non-ischemic Dilated Cardiomyopathy.

Patient-specific modelling has emerged as a tool for studying heart function, demonstrating the potential to provide non-invasive estimates of tissue passive stiffness. However, reliable use of model-derived stiffness requires sufficient model accuracy and unique estimation of model parameters. In this paper we present personalised models of cardiac mechanics, focusing on improving model accuracy, while ensuring unique parametrisation. The influence of principal model uncertainties on accuracy and parameter identifiability was systematically assessed in a group of patients with dilated cardiomyopathy ([Formula: see text]) and healthy volunteers ([Formula: see text]). For all cases, we examined three circumferentially symmetric fibre distributions and two epicardial boundary conditions. Our results demonstrated the ability of data-derived boundary conditions to improve model accuracy and highlighted the influence of the assumed fibre distribution on both model fidelity and stiffness estimates. The model personalisation pipeline-based strictly on non-invasive data-produced unique parameter estimates and satisfactory model errors for all cases, supporting the selected model assumptions. The thorough analysis performed enabled the comparison of passive parameters between volunteers and dilated cardiomyopathy patients, illustrating elevated stiffness in diseased hearts.

Analysis of lead placement optimization metrics in cardiac resynchronization therapy with computational modelling.

AIMS: The efficacy of cardiac resynchronization therapy (CRT) is known to vary considerably with pacing location, however the most effective set of metrics by which to select the optimal pacing site is not yet well understood. Computational modelling offers a powerful methodology to comprehensively test the effect of pacing location in silico and investigate how to best optimize therapy using clinically available metrics for the individual patient. METHODS AND RESULTS: Personalized computational models of cardiac electromechanics were used to perform an in silico left ventricle (LV) pacing site optimization study as part of biventricular CRT in three patient cases. Maps of response to therapy according to changes in total activation time (ΔTAT) and acute haemodynamic response (AHR) were generated and compared with preclinical metrics of electrical function, strain, stress, and mechanical work to assess their suitability for selecting the optimal pacing site. In all three patients, response to therapy was highly sensitive to pacing location, with laterobasal locations being optimal. ΔTAT and AHR were found to be correlated (ρ < -0.80), as were AHR and the preclinical activation time at the pacing site (ρ ≥ 0.73), however pacing in the last activated site did not result in the optimal response to therapy in all cases. CONCLUSION: This computational modelling study supports pacing in laterobasal locations, optimizing pacing site by minimizing paced QRS duration and pacing in regions activated late at sinus rhythm. Results demonstrate information content is redundant using multiple preclinical metrics. Of significance, the correlation of AHR with ΔTAT indicates that minimization of QRSd is a promising metric for optimization of lead placement.

Impact of coronary bifurcation morphology on wave propagation.

The branching pattern of the coronary vasculature is a key determinant of its function and plays a crucial role in shaping the pressure and velocity wave forms measured for clinical diagnosis. However, although multiple scaling laws have been proposed to characterize the branching pattern, the implications they have on wave propagation remain unassessed to date. To bridge this gap, we have developed a new theoretical framework by combining the mathematical formulation of scaling laws with the wave propagation theory in the pulsatile flow regime. This framework was then validated in multiple species using high-resolution cryomicrotome images of porcine, canine, and human coronary networks. Results demonstrate that the forward well-matchedness (no reflection for pressure/flow waves traveling from the coronary stem toward the microcirculation) is a salient feature in the coronary vasculature, and this result remains robust under many scenarios of the underlying pulse wave speed distribution assumed in the network. This result also implies a significant damping of the backward traveling waves, especially for smaller vessels (radius, <0.3 mm). Furthermore, the theoretical prediction of increasing area ratios (ratio between the area of the mother and daughter vessels) in more symmetric bifurcations found in the distal circulation was confirmed by experimental measurements. No differences were observed by clustering the vessel segments in terms of transmurality (from epicardium to endocardium) or perfusion territories (left anterior descending, left circumflex, and right coronary artery).